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. Introduction & Summary of Approach and Results 


Flare forecasting is important for protecting astronauts and space assets. 
How can forecasting be improved? 
Flares involve the largest current densities in the solar atmosphere. 


The magnetic field B, and hence current density J = cV x B/(47) are most 
accurately known at the photosphere. 


Large flares (M, X) occur in the neutral line regions (NLRs) of active regions 
(ARs). J is max. in NLRs. Byerticay changes sign. Max. free energy. 


Are there changes in the photospheric J in NLRs of ARs that are useful for 
forecasting M/X flares? 


Can HMI (Helioseismic & Magnetic Imager) be used to detect previously 
undetected changes in J? 


HMI: Full disk, continuous time observations of photospheric B at 1”, 12 


minute resolution. High enough to begin to resolve granulation dynamics: 
space and time scales ~ 1000 km and ~ 15 — 20 minutes. 


Existing flare models we are aware of do not compute the complete J. 
J #0< non-potential B. Flares relax B towards J = 0: Min. energy state. 


Objective: Develop a model that computes J using highest resolution data 
with continuous space and time coverage of entire ARs. 


Problem: Photospheric measurements return B(z,y,t), not B(z,y,z,t). So 
J, can be computed, but not J; or Jy. 


Solution: We combine HMI data, the V -B = 0 constraint, and a Fourier 
expansion of B in (x,y) to determine B(z, y, z,t) through second order in z, 
where z = 0 is the photosphere. 


Then the complete J(z,y,0,t) can be computed. 


e From B(z,y,z,t) compute J,A,E in each pixel for 14 ARs. 7 with M/X 


flares. 7 with B, C, or no flares. 


e Compute time series of Q(t) = 7J? for each AR NLR(t). Are there correla- 
tions between changes in Q and M/X flare times? 


Answer: Plausibly yes. 


e Compute the cumulative distribution function (CDF) N(Q) for each AR 
time series of Q. N(Q) is the number of events with heating rates > Q. 
Compare with the observed N(F) for the total energy FE released in solar 
flares. EF is the amount of magnetic energy converted into particle energy. 


Result: Like N(E), N(Q) is a scale invariant power law distribution: N(Q) = 
CarQ-°. N(Q) and N(E) have essentially the same exponent range. 


2. Magnetic Field 


Let L, and L, be the x and y dimensions of the rectangular region used to 
enclose most or all of the AR modeled. The HMI pixel side length A = 0.5”. 
The number of HMI data points covering this region is N = (N, + 1)(N,+ 1), 
where N, = L,/A,N, = L,/A, and (N,, N,) are given by the HMI datasets. For 
sufficiently small z, 
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e bym(t) are complex, and L(x, y,t) = Lo(z,y,t)+ 2li(2,y,t)/Lo where Lp and Ly; 
are real and determined by the HMI data and the V - B = 0 condition. 


e For z = 0, and given the N vectors B(z;, y;,0,t;) from the HMI data for each 
j, Eq. (1) is solved for the time series of the b,,,,(t) using a FFT. 


2.1. The V-B=0 Condition 


Define By = B(z,y,0,t). Take the divergence of Eq. (1) and set it equal to zero. 
Solving the resulting equation through order z gives 
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The right hand sides of Eqs. (2) and (3) are evaluated at z = 0. Therefore, Lo 
and L,; are completely determined by the HMI data. Gives 0B/0z at z = 0. 


It is the V-B = 0 condition plus the HMI data that determine the z depen- 
dence of the model that allows J, and J, to be computed. 


3. Current Density 
The current density J(x,y,z,t) =cV x B/(47). Through order z, 
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4. Vector Potential and Electric Field 
Assume the following expansion, valid through order z? for sufficiently small z. 
A(z, Y,%, t) = a(x, Y; t) se ai(x, Y; t)z oe a(x, Y; t)z? a3(x, Y; t)z°. (7) 


Expand B through order z? to obtain 
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e The a;(0 <i <3) are found by solving A = V x B with the Coulomb gauge 
condition V - A = 0 (ensures uniqueness), order by order in powers of z. 
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Then E = —c“'0A/0t — V® ~ —c7!0A/Ot. 


=0(0. Then 


5. Need to Remove Spurious Doppler Periods From the HMI B 


There is spurious, Doppler shift generated noise in the form of 6, 12, and 
24 hour period oscillations in the components of B for each pixel. Noise is 
due to SDO orbital motion. 


The noise is removed from the time series of HMI B for each pixel using 
an FFT based bandpass filter. 


Doppler Noise in Pixel Level Quantities: 


Figures 1-3 shows the filtered and un-filtered HMI time series of B,, B,, 
and B, for a randomly selected pixel from the NLR of NOAA AR 1166 
during a 70 hour long time series. 
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Figure 1: Comparison of the Filtered and Un-filtered B, in a pixel from NLR/AR 1166. 
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Figure 2: Comparison of the Filtered and Un-filtered B, in a pixel from NLR/AR 1166. 
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e Doppler Noise in NLR Integrated Quantities: 


Figures 4-7 show the results of integrating the filtered and un-filtered pixel 
level results for 7J? and B?/8z over the NLR at each time. 


The 70 hour long time interval includes 1 X, 2 M, and 9 C flares. 
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Figure 3: Comparison of the filtered and un-filtered, NLR integrated 7.J?. 
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Figure 4: Comparison of the filtered and un-filtered, NLR integrated 7.J?. 
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Figure 5: Comparison of the filtered and un-filtered, NLR integrated B?/87. 
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Figure 6: Comparison of the filtered and un-filtered, NLR integrated B?/8>7. 
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6. NLR Integrated Resistive Heating Rates Q(t) of Strongly Flaring (SF) ARs: 
Comparison with Flare Times 


e Plot Q(t). Superimpose the times of C, M, and X flares. NOAA database. 
e Look for changes in Q that may be correlated with flare times. Large spikes. 


e Check to what extent J is force-free (i.e. J_/J) <1). The largest heating 
events, which are found to occur in single pixels, are highly non-force-free: 
J x B=0 is a bad assumption at the photosphere. 


e Results for weakly flaring ARs (C, B, or no flares) are not shown. Q is 
~ 10 — 100 times less than for SF ARs, and there does not seem to be a 
correlation between spikes in Q and subsequent flaring. 
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Figure 7: NLR integrated Q for 2 of 7 SF ARs. 
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Figure 8: NLR integrated Q for 2 of 7 SF ARs. 
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Figure 9: NLR integrated Q for 2 of 7 SF ARs. 


20 


NLR Integrated nw (ergs-s-'—cm"') 


NOAA AR 2017 


18 es tte ae el ea 
reed fae reo 

O00O0 0 COOCOOO 
oo fo) 


-100 -50 50 


t (Hours Relative to X1.0 Flare at t=0. Flare day/time=20140329/1735) 


Figure 10: NLR integrated Q for the 7th SF AR. 


7. Scale Invariant Power Law Distributions of Q - Comparison with Flares 


e Compute the Cumulative Distribution Functions (CDFs) of the time series 
for Q. The CDF is the number N(Q) of heating events with heating rate 


2Q. 
e Above an AR dependent threshold value, the CDF for each AR is well fit 


by a scale invariant power law distribution N(Q) = AQ~*, with S constant 
over a range of several orders of magnitude in Q. 


Scale invariance means that a change in scale of Q (i.e. replacing Q by kQ) 
does not change the form of N(Q) (i.e. N(kQ) = constant x N(Q)). N(Q) is 
scale invariant over the range of Q for which S is constant. 
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For the 14 ARs analyzed it is found that 0.40 < S < 0.53, with a mean and 


standard deviation across the ARs of 0.47 and 0.045. This = little statistical 
variation in S from one AR to another. 


The CDF N(E£) for the total energy F released in solar flares is determined 
from observations to have the same form: N(E) = constant x E~°. 


EUV and SXR observations of nanoflares in the 0.7—4 MK range, and HXR 
observations of flares imply that 0.51 < a < 0.57, and 0.4 <a < 0.6. 


Observations also show that, as is found for the exponent S in N(Q), there 
is little variation of a among ARs. 


Therefore, the power law scaling of the photospheric Q is essentially identical 
to that found for coronal flares. 


Suggests the mechanisms that generate Q and coronal flares are closely 
related. 
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Figure 11: CDFs for all ARs, and all SF ARs. 
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8. Conclusions 


Flare forecasting models based on computing time dependent maps of the 
complete photospheric current density need to be developed. 


The spurious Doppler periods in the HMI magnetic field can introduce large 
errors into B and derived quantities. 


In combination with the model, HMI might be revealing previously unde- 
tected photospheric heating events on granulation space and time scales in 
NLRs of ARs. 


The largest heating events occur in NLRs that exhibit M/X flares. They 
are highly non-force-free. 


It is plausible these events are correlated with M/X flares, preceding them 
by several hours to several days. But the sample size of 14 ARs is too small 
to determine if a correlation exists. Analysis of more ARs is needed. 
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e The CDFs of Q obey a scale invariant power law distribution essentially 


identical to that of the energy released in flares. This suggests a close 
connection between the process that drives @, which is a photospheric 
phenomenon, and the process that drives flares, which are a coronal phe- 
nomenon. 
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